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ABSTRACT 

We conduct a systematic examination of the properties of models for chemically homogeneous, differentially 
rotating, main-sequence stars of mass 1-2 Mq. The models were constructed using a code based on a re- 
formulation of the self-consistent field method of computing the equilibrium stellar structure for a specified 
conservative internal rotation law. The code has recently been upgraded with the addition of new opacity, equa- 
tion of state, and energy generation routines, and a mixing-length treatment of convection in the outer layers of 
the stellar interior Relative to nonrotating stars of the same mass, these models all have reduced luminosities 
and effective temperatures, and flattened photospheric shapes (i.e., decreased polar radii) with equatorial radii 
that can be larger or smaller, depending on the degree of differential rotation. For a fixed ratio of the axial 
rotation rate to the surface equatorial rotation rate, increasingly rapid rotation generally deepens convective 
envelopes, shrinks convective cores, and can lead to the presence of a convective core (envelope) in a 1 M© 
(2 Mq) model, a feature that is absent in a nonrotating star of the same mass. The positions of differentially 
rotating models for a given mass M in the H-R diagram can be shifted in such a way as to approximate the non- 
rotating ZAMS over ranges in luminosity and effective temperature that correspond to a mass interval between 
M and about 0.7 M. We briefly note a few of the implications of these results, including (i) possible ambiguities 
arising from similarities between the properties of rotating and nonrotating models of different masses, (ii) a 
reduced radiative luminosity for a young, rapidly rotating Sun, (iii) the nuclear destruction of lithium and other 
light metallic species in the layers beneath an outer convective envelope, and (iv), the excitation of solar-hke 
oscillations and the operation of a solar-like hydromagnetic dynamo in some L5 - 2 Mq stars. 
Subject headings: stars: interiors — stars: rotation 
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L INTRODUCTION 

Rotation is a universal stellar physical attribute. This con- 
clusion is supported by an enormous body of observational 
data, accrued over a period of nearly 100 years, from which 
information about rotational speeds and periods, the depen- 
dence of these quantities on parameters such as mass and 
age (i.e., evolutionary state), and the effects of rotation on 
the shape, effective temperature, chemical homogeneity, and 
other basic properties of stars have been inferred. Despite in- 
creasing evidence that rotation can have a significant impact 
on a variety of stellar characteristics, it is generally not in- 
cluded as a component of the structural/evolutionary models 
which are the primary tools for the interpretation of observa- 
tions. This omission is, in part, a consequence of the increased 
complexity of the problem of determining the structure and 
evolution of a rapidly, differentially rotating star: a one- 
dimensional model necessarily becomes two-dimensional, the 
gravitational potential must be derived by solution of Pois- 
son's equation, and the uncertain physics of convective and 
circulatory flows, along with other rotation-dependent mecha- 
nisms that contribute to angular momentum redistribution and 
chemical mixing, needs to be addressed. 

Some progress toward the development of a straightforward 
yet robust technique for computing the internal structure of 
a rotating star has recently been made with the implementa- 
tion of a new version of the self-consistent field (SCF) method 
(Jackson, MacGregor, & Skumanich 2005, hereafter Paper I). 
In its original form (Ostriker & Mark 1968; see also Jack- 
son 1970), the SCF approach to computing a model for an 

* The National Center for Atmospheric Research (NCAR) is sponsored by the 
National Science Foundation. 



assumed conservative law of rotation consisted of two sep- 
arate steps: (i) determination of the gravitational potential 
for a given distribution of the mass density p, followed by 
(ii) solution of the equations of stellar structure with the po- 
tential from (i) to update the equilibrium distributions of p 
and other quantities. The process was initiated with a trial 
distribution for p, and the two steps were executed sequen- 
tially and iterated upon until the densities from (i) and (ii) 
agreed to within a specified tolerance. While the SCF method 
was successfully used to construct detailed models of rotating 
upper-main-sequence stars (Bodenheimer 1971), it failed to 
converge for objects less massive than about 9 Mq, behavior 
which precluded its application to intermediate- and low-mass 
stars (see, e.g.. Chambers 1976; Clement 1978, 1979; Paper 
I). 

As described in Paper I (see also §2 below), the refor- 
mulated SCF method circumvents the difficulties responsi- 
ble for nonconvergence in lower-mass, more centrally con- 
densed stars with an approach that entails the specification of 
trial functions for the pressure P, the temperature T , and the 
shape of constant-density surfaces, together with iterative ad- 
justment of both the profiles and central values of P and T . 
With these modifications, the method is capable of produc- 
ing self-consistent models for rotating main-sequence stars of 
all masses. It has been validated through detailed compar- 
isons with stellar models (both rotating and nonrotating) for 
masses > 2 Mq computed by other investigators using alter- 
native techniques (Paper I), and has been applied in an exam- 
ination of the photospheric shape of the Be star Achernar, an 
object revealed by interferometric observations as highly flat- 
tened by rapid rotation (Jackson, MacGregor, & Skumanich 
2004, and references therein). More recently, the structure 
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code based on the method has undergone considerable reno- 
vation, with the replacement of routines for the equation of 
state, nuclear energy generation rates, and opacities, and the 
addition of a mixing length treatment of convection. It is now 
equipped for use in an investigation of the effects of rotation 
on the structure and properties of stars with masses < 2 Mq. 

In the present paper, we use the formalism described above 
to conduct a detailed survey of the structural characteristics 
of differentially rotating stars having masses M in the range 
1 < M < 2 Mq. In the absence of rotation, this mass inter- 
val encompasses considerable variation in internal properties, 
with the gross structure of objects at the lower limit consist- 
ing of an inner radiative core and an outer convective enve- 
lope, changing to a convective core and radiative envelope for 
objects at the upper limit. The influence of rotation on the 
basic morphology of stellar interiors for these masses has not 
received much attention, with rotation-related effects usually 
treated as perturbations to the nonrotating structure, if at all 
(see, e.g., Thompson et al. 2003). However, observations in- 
dicate the occurrence of surface rotation speeds rapid enough 
to imply non-negligible modifications to many stellar proper- 
ties, particularly if the rotation is differential: the projected 
rotation speeds of 2 Mq main-sequence stars are typically in 
excess of 100 km s"^ and comparable values have been mea- 
sured for near zero-age main-sequence (ZAMS) 1 Mq stars in 
young clusters (Stauffer 1991; Wolff & Simon 1997; Tassoul 
2000, and references therein). 

Motivated by these considerations, we have constructed an 
extensive set of self-consistent models for 1-2 Mq ZAMS 
stars, to systematically study the dependence of stellar char- 
acteristics on the rate and degree of differential rotation. Our 
SCF models extend previous computational results for stars 
in this mass range, some of which were obtained for uni- 
formly rotating configurations (Faulkner, Roxburgh, & Strit- 
matter 1968; Sackmann 1970; Kippenhahn & Thomas 1970; 
Papaloizou & Whelan 1973; Roxburgh 2004), while others 
were obtained using either approximate methods for noncon- 
servative differential rotation (Endal & Sofia 1981; Pinson- 
neault et al. 1989; Eggenberger, Maeder, & Meynet 2005) or 
a non-SCF, finite-difference technique in the case of conser- 
vative rotation (Clement 1979). In §2 we provide a synopsis 
of the new SCF method, briefly describing its implementa- 
tion in a code for computing the structure of chemically ho- 
mogeneous, differentially rotating main sequence stars, the 
improvements and extensions to the input physics that have 
been made since Paper I, and the results of tests of code reli- 
ability through comparisons with extant models for nonrotat- 
ing ZAMS stars with masses < 2 Mq. An examination of the 
properties of models spanning a wide range of internal rota- 
tion characteristics is presented in §3, with particular attention 
paid to the behavior of such quantitites as luminosity and ef- 
fective temperature, central thermodynamic properties, the lo- 
cation and extent of convective regions, and the size and shape 
of the stellar photosphere. Among the consequences of these 
rotation-induced changes in stellar properties is a shift in the 
position of objects in the classical H-R diagram (HRD); in- 
deed, for some adopted rotation laws, the resulting structural 
modifications can enable a differentially rotating star to have 
the same 7^// and L values as a non-rotating star of signifi- 
cantly lower mass and to thereby occupy the same position in 
the HRD. In the concluding section of the paper (§4), we use 
the model results to consider how the changes brought on by 
rapid differential rotation might affect other features of main- 
sequence stars in this mass range, such as the radiative lumi- 



nosity of the young Sun and associated effects on the planets, 
the abundance of lithium in the surface layers of mid-F dwarf 
stars, and the presence of strong, large-scale magnetic fields 
in some A stars. 

2. THE NEW SCF METHOD 

2.1. Implementation 

The new version of the SCF method (described in detail in 
Paper I) is an iterative scheme that is initialized by specify- 
ing (i) a pair of one-dimensional trial functions, one for the 
temperature distribution and one for the pressure distribution 
(each normalized by its central value and defined over a spa- 
tial range that is normalized by the equatorial radius of the 
star), and (ii) a two-dimensional normalized function describ- 
ing the shape of the equidensity surfaces. The normalized 
two-dimensional trial density distribution, which is used as 
the source term in Poisson's equation for the gravitational po- 
tential, follows from the equation of state and the trial func- 
tions for P and T. For a conservative law of rotation, in which 
the angular velocity depends only on the perpendicular dis- 
tance from the axis of rotation, it is possible to define an effec- 
tive potential from which both the gravitational and centrifu- 
gal forces can be derived. The surfaces of constant effective 
potential (i.e., level surfaces) can then be identified and used 
to solve a set of ordinary differential equations analogous to 
the usual equations of stellar structure for nonrotating stars. 
This solution step yields updated temperature and pressure 
distributions, allowing the iterative cycle to be repeated and 
the process continued until convergence is achieved. When 
the input and output functions representing the normalized 
temperature and pressure profiles are in agreement, the two 
parameters corresponding to the central temperature (Tc) and 
central pressure (Pc) are adjusted by a Newton-Raphson tech- 
nique to bring them closer to the actual physical conditions 
at the center of the final equilibrium model. The entire pro- 
cedure (consisting of an SCF loop nested inside a Newton- 
Raphson loop) is repeated until an acceptable level of agree- 
ment between the input and output values of the two central 
parameters is attained. When used to construct models for 
differentially rotating main-sequence stars, this reformulation 
of the SCF method has been shown to converge for all masses 
in the range 0.6 < M < 30 Mq, and for values of the dimen- 
sionless rotational kinetic energy t (the ratio of the rotational 
kinetic energy to the absolute value of the gravitational po- 
tential energy of the configuration) as high as 0.10-0.12 for 
intermediate- and high-mass models, and up to nearly 0.26 for 
some fully convective, highly flattened, disk-like, 1 Mq mod- 
els. By comparison, the largest t values among the non-SCF 
models computed by Clement (1978, 1979) were f« 0.18 for 
30 Mq and < 0.12 for models in flie range 1.5 < M < 5Mq, 
whereas Bodenheimer (1971) obtained a 60 Mq SCF model 
with f « 0.24. 

As in Paper 1 (see also Jackson, MacGregor, & Skumanich 
2004), the internal rotation of each of the stellar models dis- 
cussed in §3 is given by an angular velocity distribution of the 
form 

nm = — , ^° ,2 , (1) 

l + {om/R,) 

where 05 = r sin 9, and Rg is the equatorial radius of the star. 
For a given model, the constants Qq and o; are prescribed 
parameters that characterize, respectively, the axial rotation 
rate and the ratio of the surface equatorial rotation rate to the 
axial rate (a measure of the degree of differential rotation). 
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i7(./0o = 1/(1 + cP')- In practice, the value of the first of these 
quantities is specified through the parameter 77 = fio/ with 
Qcr the equatorial angular velocity for which the magnitudes 
of the gravitational and centrifugal forces at Rg are equal. 

Rigorous answers to questions about the existence and 
uniqueness of solutions to complicated systems of integro- 
differential equations are generally very difficult to obtain. 
Although specifying the two rotational parameters a and rj 
for chemically homogeneous models of fixed mass and com- 
position does not always lead to a converged model, when it 
does, our experience with the code suggests that the model is 
unique. For a given set of (a,j]) values, using the SCF code to 
converge models of the same mass and chemical composition 
from two or more different trial models always seems to lead 
to the same final converged model. On the other hand, it can 
be demonstrated that specifying the parameters a and f^o does 
not, in general, lead to a unique solution. A minor drawback 
to use of the parameters (a, rj) in presenting the results is that 
several important global properties, including the luminosity 
L, the total angular momentum y, and the dimensionless rota- 
tional kinetic energy f , are not monotonic functions of r] when 
a is held fixed (see also §3.1). Despite this shortcoming, we 
feel that the apparent uniqueness of the models corresponding 
to particular {a,r}) values justifies the use of these quantities 
in presenting and discussing results. We emphasize that se- 
quences of models obtained by holding a fixed and varying 77 
should not be interpreted as any sort of evolutionary sequence. 
The significance of models characterized by the same value of 
a is that they have the same degree of differential rotation, that 
is, the same ile/ilo- The "half-width" of the rotation profile, 
f5i/2 = ^c/q^> does, however, change from model to model for 
constant a. 

2.2. Input Physics 

Since the publication of Paper I, the input physics for the 
SCF code described therein has been updated considerably, 
with the installation of software components that at various 
times were parts of the stellar-evolution code developed by 
Don VandenBerg at the University of Victoria. All of the 
models presented in this paper were computed for the fol- 
lowing abundances by weight of H, He, and heavy elements: 
X = 0.71 12, Y = 0.27, and Z = 0.0188. The opacities were ob- 
tained, as in VandenBerg et al. (2000), from tables of OPAL 
opacities calculated by Rogers & Iglesias (1992) and from ta- 
bles of low-temperature opacities calculated by Alexander & 
Ferguson (1994), using interpolation subroutines written by 
VandenBerg (1983). Other subroutines written by Vanden- 
Berg (1992) were utilized for the following; (i) the equation 
of state formulated by Eggleton, Faulkner, & Flannery (1973, 
EFF); (ii) nuclear energy generation rates for hydrogen burn- 
ing from Caughlan & Fowler (1988), including the effect of 
electron screening as treated by Graboske et al. (1973) for the 
case of equilibrium abundances of CNO isotopes; and (iii) a 
standard mixing-length treatment of surface convective zones 
(see, e.g.. Baker & Temesvary 1966; Kippenhahn, Weigert, & 
Hofmeister 1967). We have made no attempt to incorporate 
any of the direct effects of rotation into the adopted convec- 
tion model; instead, we have simply modified the nonrotat- 
ing mixing-length description of convection by replacing the 
local gravitational acceleration, g, with the effective gravita- 
tional acceleration, g<,// = g-n((!5)^(!5 e^s (i.e., g as reduced 
by the local centrifugal acceleration), averaged over equipo- 
tential surfaces. For all of the models, a value of 1 .9 for the 
ratio of the mixing length to the pressure scale height has been 



adopted. 

2.3. Validation 

In view of the revisions and updates that have been made 
to the SCF code of Paper I, it seems worthwhile to make 
a careful comparison of our nonrotating models for stars on 
the lower main sequence with models for the same mass ob- 
tained from a standard (nonrotating) stellar evolution code 
For this purpose, we have used the current version of the 
evolutionary code of Christensen-Dalsgaard (1982, hereafter 
referred to as the JCD code) to generate two evenly spaced 
sequences of seven models each, spanning the mass range 
0.8 < M < 2.0 Mq along the ZAMS. These models have 
chemical composition X = 0.711, Z = 0.019, quite close to 
the abundances adopted for our SCF models, and were com- 
puted using the same value (1.9) of the mixing length pa- 
rameter. The basic EFF equation of state was used to con- 
struct one of the sequences, while the other was derived with 
the so-called CEFF equation of state, a modification of the 
EFF treatment that includes Coulomb corrections (see, e.g., 
Christensen-Dalsgaard & Dappen 1992). Aside from the in- 
clusion or omission of Coulomb effects in the equation of 
state, the package of input-physics routines used to generate 
the JCD models is very similar to that installed in the SCF 
code (see, e.g., Christensen-Dalsgaard, Proffitt, & Thompson 
1993; Di Mauro & Christensen-Dalsgaard 2001). In addition 
to the intrinsic distinctions between the mathematical tech- 
niques used to obtain the two sets of models, the current ver- 
sion of the SCF code differs from the JCD code in the follow- 
ing ways: (i) the equation of state does not include treatment 
of Coulomb corrections; (ii) the photospheric pressure in the 
SCF code follows from the application of a different, sim- 
plified boundary condition (see Paper I); and (in), there is a 
sHght difference in the vintage of the nuclear energy genera- 
tion data. 

We have conducted a quantitative comparison of some of 
the important properties (R, L, Tejj, Pc, Tc, pc) of SCF, 
JCD/EFF, and JCD/CEFF models for nonrotating ZAMS stars 
of the same mass. A theoretical HRD indicating the positions 
of the various models is displayed in Figure 1. The EFF and 
CEFF main sequences (the dashed and dotted lines, respec- 
tively) are very nearly coincident, with the former models 
shifted shifted along the locus relative to the latter models, 
toward somewhat lower luminosities and effective tempera- 
tures. On the basis of the preceding discussion, we expect the 
SCF models to agree better with EFF models than with CEFF 
models, and inspection of the results plotted in Figure 1 in- 
dicate that this is generally the case. For the luminosity, the 
most sensitive of the stellar properties, the largest discrepan- 
cies occur for the 0.8 M© models, with the SCF L about 2% 
lower than that of the EFF model of the same mass; this latter 
value is, in turn, about 7% lower than L for the corresponding 
CEFF model. The luminosities of the SCF and EFF mod- 
els are essentially equal for 2 Mq, while a comparison of the 
EFF and CEFF models for that mass reveals that the CEFF 
model is about 3% more luminous. The effective tempera- 
tures and radii of the SCF models deviate from those of the 
corresponding EFF models by less than 1%, except for those 
models between 1.2 and 1.4 Mq where the the magnitudes of 
the discrepancies are somewhat larger, « l%-2%. The SCF 
values for Pc and pc differ from the EFF values of those quan- 
tities by about 2%, with the relative difference in the T^. values 
< 1%. There is also good qualitative agreement between the 
SCF and JCD models with respect to the appearance or ab- 
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Fig. 1. — A theoretical HR diagram showing the positions of nonrotating, 
chemically homogeneous stellar models for die indicated masses, as com- 
puted using different codes and equations of state. The various symbols de- 
note models obtained using the SCF code of the present paper (*), a current 
version of the code of Christensen-Dalsgaard (1982, JCD) with the stmpU- 
fied equation of state of Eggleton, Faulkner, & Flannery (1973, EFF) (A), 
and the JCD code with a version of the EFF equation of state that includes 
Coulomb corrections (CEFF) (□). The dashed and dotted Unes are the com- 
puted ZAMS locations derived from the JCD/EFF and JCD/CEFF models, 
respectively. 



sence of convective cores and envelopes, and there is good 
quantitative agreement with respect to the radial extent and 
the enclosed mass of the principal radiative-convective inter- 
faces. 

3. PROPERTIES OF DIFFERENTIALLY ROTATING STELLAR 

MODELS 

3.1. Convergence Properties 

The properties of SCF models for rotating ZAMS stars with 
masses I < M <2 Mq are summarized in Figures 2-8 and 
Table 1 . To facilitate discussion of the computed stellar char- 
acteristics, we adopt the convention, established in Paper I, of 
identifying each model by its mass, M, and the two rotational 
parameters, a and rj, for the reasons discussed in §2. Regions 
of the (a, rj) parameter space in which converged 6 Mq SCF 
models can and cannot be obtained have been illustrated in 
Figure 3 of Paper I. While the (a, 77) planes for the lower- 
mass SCF models considered here closely resemble that for 
the 6 Mq models, there are some important differences. Of 
particular relevance to the present paper are the regions cor- 
responding to Region II in Paper I, regions of relatively high 
angular momentum wherein the SCF method is incapable of 
producing converged models. These forbidden zones are con- 
siderably expanded for the 1 and 2 Mq models: for the for- 
mer models, the values (a,, t],) corresponding to the lower 
tip of this region are a, = 1.39, r], « 2.4, while for the latter, 
a, = 2.83 and 77, w 6.5. Along each constant-a sequence for 
< a < at, quantities such as the axial angular velocity fio 
and the total angular momentum / are non-monotonic func- 
tions of ij, increasing to a maximum for r/ w t], and decreas- 
ing thereafter As a ^ ai, convergence in the vicinity of this 
maximum becomes significantly slower, and fails completely 
when a is large enough (i.e., > at) to place the model within 
the forbidden zone. For a, < a < 7 (the largest value of a we 
have considered), converged models can be readily obtained 
on the low-ry sides of these regions. On the high-r; sides, 
converged models can be obtained only within the ranges 




Fig. 2. — Selected properties of differentially rotating, 1 Mq , ZAMS stellar 
models. The model characteristics are shown as functions of r; for < a < 5, 
where the parameters a and 77 specify the assumed internal angular velocity 
distribution given by equation (1). The quantities depicted in the various pan- 
els include: (a) the surface equatorial rotation speed Ve; (b) the luminosity L 
in units of Lq, the luminosity of a nonrotating 1 Mq model; (c) the central 
temperature relative to the corresponding value T^q for the nonrotating 
model; (d) the average effective temperature T^f/', (e) the equatorial radius 
Re in units of Rq, the radius of the nonrotating model; (t) the maximum per- 
pendicular distance Zmm from the equatorial plane to the photosphere (solid 
curves), and the polar radius Rp (dotted curves), as fractions of R^; (g) the 
radii of the base of the convective envelope rce and the convective core i-cc 
both measured in the equatorial plane relative to Re ; and, (h) the temperature 
Tee at the base of the convective envelope. 



1.39 < a < 3 for the 1 Mq models, and 2.83 < a < 3.47 for 
the 2 Mq models. The physical, mathematical, and computa- 
tional reasons for the lack of convergence of models on either 
side of the forbidden zone are discussed in considerable de- 
tail in Paper I. In this paper, we present results for complete 
constant-a sequences of models having a < at, but confine 
our attention to just the low-r; sides of forbidden zones for 
models with a> a,. Models on the high-?7 sides of these re- 
gions have highly rotationally flattened, disk-like structures 
that will be the focus of a subsequent paper. 

3.2. 7 Mg, Models 

In Figure 2, we show how some of the characteristics of the 
1 Mq models depend on the dimensionless rotation parame- 
ter 7] for < a < 5. The relation between rj and the physical 
quantity Vg, the equatorial rotation speed at the stellar surface, 
is given in panel (a), from which it can be seen that for each 
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Note. — Quantities listed (units in parentlieses): total mass, M {Mu,); rotational parameters, a, 77; dimensionless rotational kinetic energy, 
f; total angular momentum, J {Vf^ g cm^ s"'); equatorial velocity, Ve (km s"'); luminosity, L (Lq); mean effective temperature, (K); 

mean surface gravity, g,, (cm s"^); equatorial radius, (Rq); polar radius, Rp; maximum (normal) distance from the equatorial plane to the 
surface of star, Zmax', central pressure, P, (dyn cm""); central temperature, Tc (K); central density, p, (g cm"'); mean density, p; distance in the 
equatorial plane from the center to the top of the convective core, r„, and to the bottom of the convective envelope, fee', mass enclosed by the 
upper bounding surface of the convective core, ntcc-, and by the lower bounding surface of the convective envelope, nice; and temperature at the 
bottom of the convective envelope. Tee (K). 



of the a sequences depicted, Ve increases monotonically with 
T]. Along the curves with a < in Figure 2, the plotted mod- 
els span the range from 77 = (i.e., nonrotating) to the value 
r] = l + a^ for which Q.^ = ilcr- Along the curves with a> a,, 
the last plotted model is located adjacent to the boundary of 
the forbidden region in the (a,!]) plane; for these differen- 
tially rotating models, the centrifugal and gravitational forces 
have nearly equal magnitudes in the core of the star. We note 
that the degree of differential rotation increases with a, in the 
sense that the configurations corresponding to larger values 
of a have a greater difference between the axial and surface 
equatorial rates of rotation (see §2.1). 

As panel (b) of Figure 2 makes evident, the radiative lumi- 
nosities of these models are diminished relative to the lumi- 
nosity L() of a nonrotating 1 Mq star. This is a well-known 
consequence of including rotation in the determination of the 
equilibrium stellar structure (see, e.g., Clement 1979; Boden- 
heimer 1971). In the results shown in Figure 2, the reduc- 
tion in L is larger for differentially rotating models than it 
is for models that are uniformly or nearly uniformly rotat- 
ing. A model with a = rotating at the break-up rate (77 = 1) 
has L/Lq = 0.78, while an a = 2 model with rj = 2.42 has 
L/Lo = 0.15,a reduction of more than a factor of 6 from the 
nonrotating value. Much of the reason for this behavior lies 
in the effect of rotation on the thermodynamic conditions in 
the deep, energy-producing regions of the stellar interior. For 
these 1 Mq models, the contribution of the centrifugal force to 
supporting material against gravity enables the star to emulate 
an object of lower mass with correspondingly reduced values 
of Pc, Tc, and pc (e.g., Sackmann 1970). The results presented 
in panel (c) illustrate the dependence of Tc on model rota- 



tional properties; similar variations are found for both Pc and 
pc. For rigidly rotating configurations, this centrifugal sup- 
port is largest in the outermost layers of the interior, which 
contain only a small fraction of the stellar mass; in this case, 
Pc, Tc, and pc are little changed from the values appropri- 
ate to a nonrotating star of the same mass. For the a = 0, 
rj = I model noted previously, Pc/PcO = 0.94, Tc/Tco = 0.96, 
and pc/PcO = 0.98, where the subscript indicates the nonro- 
tating value. Alternatively, in models for higher values of a, 
the effects of rotation are increasingly concentrated toward the 
central regions of the star, with the result that the perturbations 
to the central thermodynamic quantities can be more substan- 
tial; for a = 2, = 2.42, Pc/Pco = 0.56, Tc/Tco = 0.68, and 
Pc/PcO = 0.81. Panels (b) and (c) also indicate that the mag- 
nitudes of the changes in L, Tc, and other quantities depend 
on the assumed profile of internal differential rotation. The 
model for a = 5, 77 = 4 has L/Lq = 0.54, with Pc/PcO = 0.72, 
Tc/Tco = 0.84, and Pc/pco = 0.85, smaller reductions relative 
to the nonrotating model than those for a = 2, = 2.42. This 
behavior is an outgrowth of the structural modifications aris- 
ing from the centrifugal force distributions associated with the 
different rotation profiles. In the a = 5 model, the ratio of 
the centrifugal to gravitational force in the equatorial plane, 
^r/g (r is the radial coordinate in the equatorial plane), is 
sharply peaked in the innermost portion of the stellar core, 
with maximum value 0.83 at the center, decreasing to « 0.1 
atr/Re = 03. For the shallower angular velocity profile of the 
a = 2 model, the force ratio decreases from a smaller central 
value of 0.42 to 0.23 at the stellar surface, 10 times the value 
found throughout the outer 50% of the interior of the a = 5 
model. 
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In this connection, we note that Bodenheimer (1971) found 
that the luminosities of models for rotating 30 Mq stars, com- 
puted assuming a variety of internal rotation laws, depended 
primarily on the total angular momentum content of a given 
stellar model and not on the details of its distribution within 
the interior. Specifically, his results indicated that the lumi- 
nosities of models corresponding to four different prescribed 
distributions of the angular momentum per unit mass de- 
creased with increasing total angular momentum J, with the 
relation between L and J nearly the same for each model se- 
quence. If the luminosities of our 1 Mq models are plotted 
versus their respective / values, an analogous reduction in L 
for increasing J can be discerned. However, the relation be- 
tween L and J is roughly independent of a only for small J, 
and exhibits a clear dependence on a that becomes increas- 
ingly pronounced as J is made larger. That the luminosity can 
be even approximately expressed as a function of J necessar- 
ily reflects the modifications produced by rotation to the ther- 
modynamic conditions in the energy-producing core of the 
star, as described previously (see also Mark 1968). The na- 
ture and origin of the relation between L and J for stars with 
masses 1-2 Mq and higher, including its dependence on the 
structural characteristics of the models, will be addressed in 
detail in a subsequent paper in this series. 

As in Paper 1, we define an average effective temperature 

through the relation Tejj = (L/ctA) where a is the Stefan- 
Boltzmann constant and A the area of the stellar surface. The 
results for 1^// depicted in panel (d) exhibit dependences on a 
and 7] that are similar to those seen for the quantities plotted in 
the preceding panels. Reductions by more than 1500 K from 
the nonrotating value (Tgff = 5540 K) are possible as, for ex- 
ample, in the case a = 1.5, 77 = 2.18, for which T^// = 3880 A". 
The rotation-induced variations in T^ff represent the com- 
bined effects of changes in both L and A. An indication as 
to the behavior of A can be gleaned from examination of the 
influence of rotation on the stellar size and shape. This in- 
formation is presented in panels (e) and (f), where we show, 
respectively, the equatorial radius Re (in units of the radius 
Rq of the nonrotating, spherical 1 Mq model), and the polar 
radius Rp together with Z„„,f , the maximum perpendicular dis- 
tance from the equatorial plane to the stellar surface. These 
latter two quantities are both given as fractions of Re', values of 
Rp/Re that are < 1 reflect a rotational flattening of the config- 
uration, while values of Zmax/Re that are ^ Rp/Re are indica- 
tive of a deviation from a convex spheroidal shape through the 
development of a concavity in each of the two polar regions 
of the star 

Note that the response of the equatorial radius to increasing 
77 differs greatly depending upon whether a is < 1 .5 or > 2. In 
the former case, the centrifugal force attains its largest value 
relative to the gravitational force at the stellar surface, conse- 
quently producing a distension of the outer, equatorial layers 
of the stellar interior and an overall increase in Although 
Rp decreases somewhat relative to the radius of the corre- 
sponding nonrotating model, the net effect of the changes in 
Re and Rp is usually an increase in the volume of the rotating 
star. For some models with a w a,, the decrease in Rp can 
more than compensate for the the increase in Re, and the vol- 
ume of the rotating star is reduced. In the case where a > 2, 
both Re and the steUar volume shrink with increasing t] along 
a constant-a sequence. For these models, the centrifugal force 
is largest in comparison to the gravitational force in the central 
region of the stellar core, causing the central thermodynamic 



conditions, luminosity, and size (i.e, radius, surface area, vol- 
ume) to assume values that are characteristic of nonrotating 
stars of lower mass. This raises the possibility that a differen- 
tially rotating star can imitate a less massive nonrotating star 
in radius and effective temperature, as well as in luminosity. 
In panel (f), it can be seen that for a < 1.5, the curves for 
'Z'max/Re (solid liucs) and Rp/Re (dotted lines) are coincident 
and < 1 , implying that the photospheric shape of these models 
is oblate spheroidal. Such a model (a = 1.5, 77 = 1.55) is de- 
picted in panel (a) of Figure 3. For a > 2, however, Z,^ > Rp, 
symptomatic of the development of a "dimple" or indentation 
at either pole, as in the cases of the models shown in panels 
(b) (a = 3.75, r] = 3.58) and (c) (a = 5,r) = 4.15) of Figure 3. 

Panels (g) and (h) of Figure 2 contain results pertaining to 
the location, extent, and properties of convective regions in 
the models. In the absence of rotation, the internal structure 
of a 1 Mq star consists of an inner, radiative core that encom- 
passes « 72% of the stellar radius and contains « 97% of the 
stellar mass, surrounded by an outer, convective envelope. As 
is apparent in panel (g), the inclusion of rotational effects can 
modify this basic morphological picture in two ways: either 
by increasing the size of the convective envelope or by pro- 
moting the formation of a convective core. 

For lower values of a (i.e., a < 3 in panel [g]), the rota- 
tionally induced changes in internal structure cause the outer 
convection zone to deepen as r] increases. For example, in 
the model with a= 1.5, rj = 1.55 (see panel [a] in Figure 3), 
the radius in the equatorial plane of the base of the convective 
envelope is rce/Re = 0.618, significantly deeper than the the 
base radius rce/Re = 0.718 in the nonrotating model. Associ- 
ated with this reduction in is a decrease in the mass of the 
core (nic/M « 0.90), and an increase in the temperature Tee 
at the bottom of the envelope, as can been seen in panel (h); 
for the a = 1.5, 7/ = 1.55 model, T^e = 2.93 x lO*" K, as com- 
pared with Tee = 2.46 X 10^ A" for the nonrotating model. Since 
some chemical species (e.g., Li, Be, and B) can be destroyed 
by thermonuclear reactions at temperatures > 2.5 x 10*" ^, en- 
hancements of Tee of this magnitude are likely to have conse- 
quences for the surface abundances of these elements. In the 
case of the strongly differentially rotating models with a > 3, 
the largest centrifugal effects are concentrated in the inner- 
most portion of the core, so that the fractional thickness of 
the outer envelope is little affected. However, as the magni- 
tude of the centrifugal-to-gravitational force ratio in the core 
increases for larger 77, the decreasing pressure gradient im- 
pUed by the requirement of hydrostatic equilibrium forces the 
otherwise stably stratified central regions of the interior to be- 
come convective. Under these conditions, with the presence 
of a convective core, the structure of the deep interior resem- 
bles that of a higher-mass star. The models shown in panels 
(b) and (c) of Figure 3 are examples of 1 Mq stars with con- 
vective cores; in each of these models, the radial extent of this 
region in the equatorial plane is about 14% of Re and contains 
« 2-3% of the stellar mass. For comparison, the convective 
core in a model for a nonrotating 2 Mq star has a radius that 
is about 12% of Re and contains « 14% of the steUar mass. 

To illustrate the structural differences between differen- 
tially and near-uniformly rotating stars, in Figure 4 we show 
the profiles of several physical quantities for models with 
(a, 77) = (3.00,3.26), and (1.00,1.20), along with the corre- 
sponding results for the nonrotating model. The luminosi- 
ties and total angular momenta of the two rotating models 
are L/Lq = 0.201 and 7 = 9.23 (a = 3), and L/Lq = 0.634 and 
/ = 5.30 (a = 1), respectively, where / is meaasured in units 
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(a) IMo O.93R0 




(b) IM0 O.74R0 




(c) IM0 O.77R0 




(d) I.2M0 O.9OR0 



(e) 2M0 2.93R0 





(f) 2M0 I.5IR0 




X 



Fig. 3. — Contours of level surfaces in the meridional plane for some of the nonspherical models listed in Table 1. The six rotating models shown are defined 
by the total mass and the two rotational parameters (M, a, r)) as follows: (a) 1 Mq, 1.5, 1.55; (b) 1 Mq, 3.75, 3.58; (c) 1 Mo, 5, 4.15; (d) 1.2 M©, 4, 3.74; (e) 
2 Mq, 3, 5.64; and, (f) 2 M , 4.75, 5.9. From the surface inward, the level surfaces depicted in each panel enclose a fraction of the total mass equal to 1.000, 
0.995, 0.950, and 0.500, respectively. The fractional radii in the equatorial plane of these level surfaces for the various models are: (a) 1.00, 0.88, 0.71, 0.31; (b) 
1.00, 0.90, 0.75, 0.40; (c) 1.00, 0.89, 0.72, 0.37; (d) 1.00, 0.87, 0.68, 0.34; (e) 1.00, 0.71, 0.37, 0.12; and, (f) 1.00, 0.66, 0.47, 0.24. Radiative portions of the 
interior are indicated in white, and convective regions are shaded gray. The fractional equatorial radii and enclosed masses for the interfaces between radiative 
and convective zones in the models are listed in Table 1. The nimibers at the tops of the panels denote the total mass M and equatorial radius Re of each model. 



of 10^" g cm^ s"'. In panels (a) and (b), the temperature T 
and mass density p are depicted as functions of the radial po- 
sition r (measured in units of the present-day solar radius Rq) 
in the equatorial plane of the star. The central values of both 
quantities exhibit rotation-induced reductions relative to the 
nonrotating case, the magnitudes of these modifications being 
larger for a = 3 (Tc/Td) = 0.68, Pc/pd) = 0.76) than for a = I 
(Tc/Tco = 0.92, Pc/pci) = 0.97). As noted previously, this be- 
havior is a consequence of differences in the magnitudes and 
distributions of the centrifugal force in the two rotating mod- 
els. These distinctions can be clearly seen in the profiles of the 
centrifugal-to-gravitational force ratio shown in panel (d). For 
a = 3, the force ratio is largest in the deep interior, attaining 
a maximum value of 0.96 at the stellar center and decreas- 
ing outwards to a magnitude w 0.1 in the photosphere. For 
a = 1, the ratio increases monotonically throughout the inte- 
rior, rising from a central value of just 0.026 to a maximum 
of 0.36 at Rg. The substantial contribution of the centrifu- 
gal force to the support of the innermost regions of the a = 3 
model is responsible for the considerable enhancement of the 
density scale height there, evident in panel (b); p declines by 
only ?a 10% over the iimer 20% of the stellar interior. The 
resulting changes in the internal mass distribution (panel [c]) 



lead to a star with a smaller radius, R^ = 0.72 Rq as opposed 
to Re = 0.88 Rq for the nonrotating model. Alternatively, for 
the a = 1 model, the lack of centrifugal support in the core 
of the star leads to temperature, density and mass distribu- 
tions therein that closely resemble those of the nonrotating 
model. Closer to the surface, however, the density distribu- 
tion becomes extended, a product of the increasing centrifu- 
gal reduction of gravity in the outer layers of the interior; as 
a result, the stellar radius is larger than that of the nonrotating 
model, = 1.00 Rq. 

3.3. Solar Look-Alike Models 

Figure 5 is a theoretical HRD for 1 Mq models that span a 
broad range of internal rotational states, from uniform rotation 
to extreme differential rotation (up to Qq = 50 for a = 7). 
The nonrotating ZAMS is delineated by a dotted line, with the 
positions of several specific models for 0.6 <M < 1 .0 Mq in- 
dicated. Clearly, as already implied by panels (b) and (c) of 
Figure 2, the locations of rotating models in such a diagram 
are displaced to the right of and below the position they would 
occupy in the absence of rotation, toward lower values of both 
the luminosity and the effective temperature. Models with 
uniform or near-uniform internal rotation (i.e., a = 0, 1) lie 
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Fig. 4. — Profiles of selected physical quantities in the interiors of 1 Mq 
models with ia,ri) = (3.00,3.26) (solid curves), (1.00, 1.20) (dashed curves), 
and (0.00,0.00) (dotted curves). The profiles depict the dependence of each 
quantity on radial position r in the equatorial plane, from the center of the 
star to the surface. The various panels show (a) the temperature T, (b) the 
mass density p, (c) the fraction m/M of the total mass contained within a 
level siuface that intersects the equatorial plane at r, and (d) the ratio (Cl^r/g) 
of the centrifugal force to the gravitational force. 



well above the nonrotating ZAMS, while those for which the 
degree of differential rotation is substantial {a > 3) have po- 
sitions in the HRD that collectively approximate the nonrotat- 
ing ZAMS over virtually the entire mass range depicted. For 
example, the rotating 1 Mq model (a = 3.75, r] = 3.58) shown 
in panel (b) of Figure 3 is characterized by Tejj = 4700 K 
and L/Lq = 0.224, which is indistinguishable from the val- 
ues T,jj = 4710 K, L/Lq = 0.224 for a nonrotating, 0.8 Mq 
model. Such a coincidence of the positions of stars of differ- 
ent mass in the classical HRD is made possible by the fact 
that the largest structural changes take place in the cores of 
strongly differentially rotating models. Specifically, the de- 
crease in L brought about by the smaller values of Pc, Tc, and 
pc, together with the corresponding reduction in Rg, enable 
such a model to effectively mimic a lower mass stellar model 
in which rotational effects are not included. 

Figure 6 illustrates some of the internal properties of three 
so-called solar look-alike models, that is, models for differ- 
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Fig. 5. — A theoretical HR diagram showing the positions of models for 
rotating, ZAMS, 1 Mq stars. Luminosities are given in units of the luminos- 
ity Lfl = 0.653 L0 of a non-rotating, 1 Mq star. The various symbols denote 
models constructed using the rotation law of equation (1) with the values of 
the parameters (a, r?) listed in the Figure. The ZAMS for non-rotating stars is 
indicated by the dotted line, with the positions of models for masses 0.6, 0.7, 
0.8, 0.9, and 1.0 along it marked by an * symbol. 



entially rotating stars with masses > 1 Mq that have many 
physical attributes in common with the model for a nonrotat- 
ing, 1 Mq star Inspection of panels (a) and (b) reveals that 
the profiles and the central values of the temperature and den- 
sity for the 1.1, 1.2, and 1.3 Mq models depicted therein are 
quite close to those of the solar-mass model with a = r] = ^. 
A consequence of these structural similarities is that many 
of the general characteristics of the models are nearly iden- 
tical. For example, the 1.2 M(q model which, in the absence 
of rotation, would have L/L© = 1.591, RejRQ = 1.129, and 
T,ff = 6110/: instead has L/Lq = 0.653, Re/RQ = 0.904, and 
T,ff = 5546 K, compared with L/Lq = 0.653, Re/RQ = 0.878, 
and Tgff = 5545 K for the nonrotating, 1 Mq model. A cross- 
section in the meridional plane of the rotating 1.2 Mq inte- 
rior is shown in panel (d) of Figure 3, from which it can be 
seen that a solar-like convection zone having r^e/Re = 0.746 
is present in the layers beneath the photosphere. With a sur- 
face equatorial rotation speed = 112 km s~\ such an object 
could be mistakenly identified as a rapidly rotating 1 Mq star, 
if observations were analyzed through comparison with non- 
rotating stellar models. That the values of measureable or 
inferrable stellar properties can be the same in rotating and 
nonrotating models for different masses represents a potential 
source of ambiguity in the interpretation of a variety of obser- 
vations. Panels (c) and (d) reveal subtle differences in the 
internal mass distributions and in the variation of the ratios of 
the centrifugal to gravitational force in the equatorial planes 
of these solar look-alike models. Future space-based astero- 
seismological observations may be capable of exploiting such 
differences in internal structure to distinguish slowly rotating 
stars from more rapidly, differentially rotating higher-mass 
stars that happen to have the same values of L and T^ff and 
thus occupy the same position in the HRD (see, e.g., Lochard 
etal. 2005). 

3.4. 2 Mq Models 

In Figure 7, we present a summary of the rotational depen- 
dence of 2 Mq model characteristics, using the same format as 
adopted for Figure 2. The model for a non-rotating star of this 
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Fig. 6. — Profiles in tlie equatorial plane of the temperature (a), den- 
sity (b), mass fraction (c), and the ratio of the centrifugal to gravitational 
acceleration (d) for a nonrotating, 1 Mq model (a = rj = 0) and three so- 
lar look-alike models. The latter models were obtained for (M,a,rj) = 
(1.1 M0, 4.00, 2.74), (1.2 M0, 4.00,3.74) and(1.3 M0, 4.00,4.39), andhave 
values of T^ff and L that differ from those of the non-rotating, solar model 
by about 1% or less. 



mass has Lq = 15.47 Lq,Rq = 1.59 Rq, and Teff = 9090 K, and 
a convective core with = . 1 22, nice /M = 0.141; such a 

model has a thin, subsurface convective layer, rce/Re = 0.99, 
that contains a negligible fraction of the stellar mass. The re- 
sults shown in Figure 7, which span the rotational parameter 
ranges 0<Q;<5,0<r7<7, exhibit many of the same behav- 
iors as noted previously in connection with the 1 Mq models 
of Figure 2. In particular, along each constant-a sequence, the 
luminosity (panel [b]), central temperature, (panel |c|) and ef- 
fective temperature (panel [d]) all initially decrease as rj is in- 
creased, with the largest reductions in these quantities occur- 
ring in the cases of differentially rotating models with a > 3. 
For a= 1,2, both L/Lo and Tc/Tcq pass through minima and 
then increase slightly with rj beyond those points, whereas 
Tgff continues to decrease with t] for all a sequences. The 
central pressure (panel [e]), on the other hand, is an increas- 
ing function of rj for < a < 3, and displays a non-montonic 
dependence on 77 for the model sequences corresponding to 
a = 4, 5. The origin of this variation is the mass dependence 
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Fig. 7. — Selected properties of differentially rotating, IMq.ZAMS stellar 
models, as in Figure 2. The quantities depicted in the various panels include: 
(a) Ve, (b) L, (c) Tc, (d) T,ff, (e) Pc, (t) Re, (g) Z™„ and Rp, and (h) r^e and 

fee- 



of Pc for stars on the non-rotating ZAMS. In our SCF models 
without rotation, Pc has a maximum value for M near 1 .6 Mq, 
and is decreasing for both larger and smaller M values. The 
maximum occurs for a mass close to the value marking the in- 
ternal structural transition between stars with radiative cores 
and convective envelopes and stars with convective cores and 
radiative envelopes. Hence, insofar as the properties of a ro- 
tating star are like those of a non-rotating star of lower mass, 
we expect the Pc values of these 2 M© models to increase as 
the effects of rotation become more pronounced. For cases in 
which the central thermodynamic conditions are significantly 
perturbed by rapid, differential rotation in the core of the star, 
this mass-lowering effect can be large enough to shift Pc to 
values characteristic of nonrotating stars with M < 1 .6 Mq 
(i.e., on the other side of the central pressure peak), as in the 
2 Mq models for a = 4, 5. 

Panels (f) and (g) of Figure 7 convey information pertain- 
ing to the photospheric sizes and shapes of the models. Those 
that rotate nearly rigidly are equatorially distended, and, in 
the direction perpendicular to the equatorial plane, have their 
largest dimension along the rotation axis (Z^ax = Pp)', since 
Rp < Re, these models (like their 1 Mq counterparts) have a 
flattened, spheroidal shape. Models with increasing degrees 
of differential rotation (i.e., with a > 3) develop polar con- 
cavities (as indicated by Z^ax > Rp), and ultimately (i.e., for 
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a = 4, 5 and sufficiently large 77) become more compact, with 
Re < Ro- Panels (e) and (f) of Figure 3 give cross-sectional 
representations of the photospheric shapes of 2 Mq mod- 
els for (a, rj) = (3.00, 5.64) and 4.75, 5.90), repectively; in 
the former model, Rg = 2.93 Rq > Rq, while for the latter, 

Re=l.51RQ<RQ. 

The influence of rotation on the occurrence and extent of 
convective regions within the models is explored in panel (h) 
of Figure 7. The convective core, a salient feature of the 
non-rotating stellar interior, generally decreases in size as the 
value of 7] is increased. This contraction stems from a re- 
duction in the magnitude of the radiative gradient Wrad in the 
core region of the star, a result of the way in which the cen- 
tral thermodynamic conditions are modified by rotation (see, 
e.g., MacGregor & GillUand 1986); with Vrad smaller, the 
size of the region wherein it exceeds the adiabatic gradient 
Vcui shrinks accordingly. Only for the most rapidly, differen- 
tially rotating models does I'cc show a modest increase com- 
pared to the non-rotating model; in panel [h], the model for 
(a, 7]) = (5.00,6.88) has rcc/Re = 0.124. In this case, as was 
seen for the 1 Mq models having convective cores in Fig- 
ure 2, the growth of r^c can be traced to a larger value of 
Vrarf, produced by a centrifugal- to-gravitational force ratio 
that is nearly unity at the stellar center. Note also that for 
some 2 Mq models, the convection zone underlying the stel- 
lar photosphere can extend into the stellar interior by more 
than the few 10~^ Ro that is the thickness of this region in the 
absence of rotation. As can be seen, for example, in panel 
(e) of Figure 3, the base of the convective envelope in the 
equatorial plane of the model with (a, rj) = (3.00,5.64) is lo- 
cated at rce/Re = 0.718, a fractional depth which is the same 
as in the non-rotating 1 Mq model. This leads to the intrigu- 
ing possibility that solar-like oscillations, driven by turbulent 
convection, may be excited in stars that would normally be too 
massive to generate them. The results shown in panel (h) sug- 
gest that such a solar-like convective envelope is most likely to 
occur in 2 Mq models with intermediate differential rotation 
(say., aw at= 2.83). For a profile of this kind, the inner and 
outer portions of the interior can each rotate rapidly enough 
to both significantly perturb the thermodynamic conditions in 
the core and increase Re by extending the stellar envelope. 

In Figure 8, we show the positions of models for differen- 
tially rotating, 2 Mq stars in a theoretical HRD. As in Figure 
5, the location of the non-rotating ZAMS is indicated by a dot- 
ted Une, with the positions of models for stars with masses in 
the range 1 .0 < M < 2.0 Mq indicated along it. This HRD for 
rotating 2 Mq models has a number of features in common 
with the corresponding representation of 1 Mq properties. 
Without exception, models are shifted from the non-rotating 
location to new positions characterized by lower Teff and L 
values. For a < 3, these new positions lie to the right-hand 
side of (i.e., above) the non-rotating ZAMS; for rotation that 
is increasingly differential, however, the model locations ap- 
proach the non-rotating ZAMS, moving just to the left-hand 
side of (i.e., below) it for a = 5. For a = 4.75, the model 
positions are distributed along the non-rotating ZAMS, de- 
lineating its track in the HRD for masses greater than about 
1.5 Mq. This coincidence again raises the possibility that 
a rapidly, differentially rotating 2 Mq star could effectively 
masquerade as a non-rotating star of lower mass. As an ex- 
ample, with L/Lq = 6.083, Re/Re = 1 .5 12, and Teff = 7480 K 
the model for {a, rj) = (4.75,5.90) shown in panel (f) of Fig- 
ure 3 closely resembles a non-rotating 1 .6 Mq star, for which 
L/Lq = 6.075, Re/Re = 1.466, and Teff = 7490 K. 
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Fig. 8. — A theoretical HR diagram showing the positions of models for 
rotating, 2 Mq ZAMS stars, as in Figure 5. Luminosities are given in units of 
the luminosity Lo = 15.468 Lq of a nonrotating, 2 Mq star. The ZAMS for 
nonrotating stars is indicated by the dotted line, with the positions of models 
for masses 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.8, and 2.0 Mq along it marked 
by an * symbol. 



4. SUMMARY AND DISCUSSION 

The SCF approach to computing the structure of rotating 
stars is robust and efficient, capable of yielding converged 
models of stars of all masses, with internal angular-velocity 
distributions covering the range from uniform to extreme dif- 
ferential rotation. The code described in Paper I has been up- 
graded throught the introduction of opacities, energy gener- 
ation rates, and an equation of state that are close to state- 
of-the-art, and with the introduction of a mixing-length treat- 
ment of convective energy transport. We have used this re- 
vised code to conduct a systematic study of the properties of 
models for rotating ZAMS stars with masses 1 < M < 2 Mq, 
assuming that the rotation of the stellar interior can be ad- 
equately represented by the parameterized angular-velocity 
profile given in equation (1). 

Our results suggest that rotation affects virtually every char- 
acteristic of the intermediate- and low-mass stars considered 
herein, with the magnitude, and sometimes even the sense, 
of the changes in fundamental attributes depending on both 
the rate and degree of the differential rotation in the stellar 
interior. For models that are uniformly or nearly uniformly 
rotating, the centrifugal contribution to the hydrostatic sup- 
port of the star is largest in the outer layers of the interior; 
for models that are strongly differentially rotating, the cen- 
trifugal force is largest compared to gravity in the core of the 
star. Luminosities and effective temperatures are always di- 
minished relative to nonrotating stars of the same mass, while 
equatorial radii increase in uniformly rotating or moderately 
differentially rotating models, and decrease in models where 
the difference between the angular velocity at the center and 
that at the equator is sufficiently large. As the degree of dif- 
ferential rotation is increased, the photospheric shape of the 
star changes from a convex surface closely approximating a 
spheroid flattened along the rotation axis to a roughly oblate 
surface with deepening polar concavities; in such cases, the 
greatest perpendicular distance from the equatorial plane to 
the surface occurs away from the axis of rotation, and is al- 
ways less than the equatorial radius of the star. In 1 Mq mod- 
els, the effects of rotation can either increase the thickness of 
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the outer convective envelope, or contribute to the formation 
of a convective core. In rotating 2 Mq models, the size of the 
convective core is diminished relative to that found in the ab- 
sence of rotation, and an extensive, solar-like convection zone 
can be present in the outer layers of what would otherwise be 
a stable radiative envelope. 

The extent to which the results of these computations accu- 
rately depict the stellar structural modifications arising from 
rotation depends on the validity of a number of simplifying 
assumptions and approximations made in the course of devel- 
oping the basic model, as well as on the inherent limitations 
of the SCF method itself. Here, we briefly address the most 
salient of the factors that might restrict the applicabihty of 
these results, noting where improvements and extensions can 
(or cannot) be made (see also Paper I). 

1) . The modified SCF iterative scheme utilized in the 
present investigation necessarily requires the angular velocity 
distribution within the star to be conservative, and thus ex- 
pressible as a function of just the perpendicular distance from 
the rotation axis, Q = fi(OJ) (see eq. [1]). Although this speci- 
fication facilitates the construction of models through the con- 
siderable mathematical simplification it introduces, its effect 
on some computed stellar structural characteristics, such as 
the photospheric shapes and core physical properties of very 
rapidly, differentially rotating models, is likely to differ from 
that produced by a rotation law of (say) the form fl = ri(r). 
There is some evidence, from both simulations and observa- 
tions, for the occurrence of differential rotation of the type 
given by equation (1) (see, e.g., Dobler, Stix, & Brandenburg 
2006, and references therein). Yet analyses of solar acous- 
tic oscillations yield a picture of the large-scale internal dy- 
namics of the Sun which is not in accord with the angular 
velocity being constant on cylindrical surfaces (see Thomp- 
son et al. 2003). We note that use of the SCF method pre- 
cludes consideration of non-conservative rotation laws (e.g., 
fl = n(r)); however, detailed treatment of the hydrodynamical 
and magnetohydrodynamical processes affecting the internal 
rotational states of the Sun and stars is presently beyond the 
scope of any of the extant structural/evolutionary models. 

2) . The model treats only the prescribed rotational motion 
of the stellar interior, omitting any meridional circulatory flow 
and its consequent effects on the internal angular momentum 
distribution. For models with radiative envelopes, rotation- 
ally driven circulation will cause the rotation profile to devi- 
ate over time from the state given by equation (1), unless such 
evolution is mitigated by additional angular momentum trans- 
port mechansims (see, e.g., Maeder & Meynet 2000, and ref- 
erences therein). For models with convective envelopes, the 
rotation profile in the outer layers of the interior is the prod- 
uct of the complex interplay between meridional circulation 
and turbulent heat and angular momentum transport (Rempel 
2005; Miesch, Brun, & Toomre 2006); whether fl is solar- 
like or constant on cylindrical surfaces depends on the lati- 
tudinal entropy distribution in the subadiabaticaUy stratified 
layers below the convection zone. 

3) . The model employs a simplified treatment of con- 
vection, locating unstable regions by application of the 
Schwarzschild criterion and utilizing an averaged, rotationally 
modified mixing-length description (see §2.2) to determine 
the structure of an outer convective envelope, if present. Use 
of the Solberg-H0iland condition (e.g., Ledoux 1965; Kip- 
penhahn & Weigert 1990) to ascertain the onset of convective 
instability would account for the direct influence of the cen- 
trifugal force. Since models computed for the rotation law of 



equation (1) have dj/d^ > 0, where j = 5103^ is the specific 
angular momentum (see Paper I), this change would likely de- 
crease the equatorial-plane thickness of an outer convection 
zone, and produce a latitudinal variation in rce- However, for 
rapid rotation, there is considerable uncertainty regardless of 
the convection criterion/model adopted, because of both the 
rudimentary nature of extant treatments of convective energy 
transport and the universal assumption of axisymmetry among 
rotating stellar models. 

We reiterate that the present model describes a chemi- 
cally homogeneous, ZAMS star, and neglects effects associ- 
ated with the structural and compositional evolution of the 
stellar interior. We are presently developing a mean-field 
hydrodynamics-based treatment of turbulent chemical and an- 
gular momentum transport which maintains conservative ro- 
tation profiles, thus making it possible to investigate the main- 
sequence evolution of these SCF models. Despite the short- 
comings enumerated above, we believe that the basic model 
and method described herein compare favorably with alterna- 
tive approaches to determining the structure and evolution of 
rotating stars. The most widely used of these (e.g., Meynet 
& Maeder 1997) relies on the use of approximate. Roche- 
like equipotentials to represent the internal gravity of the star, 
thereby strictly Umiting it to describing slowly rotating stars 
to ensure the accuracy of the computed models. Whatever 
its drawbacks, the SCF method yields two-dimensional, ax- 
isymmetric configurations that represent fully consistent so- 
lutions to both the set of stellar structure equations and Pois- 
son's equation for the gravitational potential. 

The results presented herein have a number of imphcations 
for young stars with masses between 1 and 2 Mq. Mea- 
surements of projected equatorial rotation speeds in excess of 
100 km s"' for some solar- type stars in young open clusters 
(see, e.g., Stauffer 1991) raise the possibility that the struc- 
ture and properties of these objects could be significantly al- 
tered from those of nonrotating stars of the same mass. Such 
stars, provided that their interior rotation is sufficiently differ- 
ential, could in actuality be somewhat more massive objects 
for which surface rotation as rapid as that indicated by obser- 
vations is the normaUy expected ZAMS state. The question of 
whether or not the kind of strong differential rotation required 
to produce such ambiguity is present within the interiors of 
some low- and intermediate-mass stars may ultimately be re- 
solved through space-based asteroseismological observations, 
which should allow low-resolution inversions of the rotation 
profiles in the inner ~ 30% of the stellar radius (see Gough 
& Kosovichev 1993). Asteroseismology from space may also 
afford the means for identifying individual stars whose rota- 
tion enables them to pose as lower-mass objects, since, when 
effects associated with asphericity are not too large, the aver- 
age mode frequency spacing is sensitive to the mean density, 
a quantity which Table 1 reveals to be different for look-alike 
models. 

If rapid differential rotation is a possibility for objects in the 
mass range spanned by the models of §3, then the associated 
changes in structure and properties could have consequences 
for a variety of important astrophysical processes that take 
place within and around such stars. The reduced radiative 
luminosity of a rapidly rotating young Sun would likely influ- 
ence the evolution of the solar nebula, and further exacerbate 
the discrepancies between the properties of nonrotating mod- 
els and observational inferences indicating a higher ZAMS lu- 
nunosity (see, e.g., Sackmann & Boothroyd 2003). The rota- 
tionaUy induced deepening of a sub-photospheric convection 
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zone, together with the increase in temperature of the mate- 
rial at the base, could contribute to the depletion of lithium 
in the stellar surface layers by reducing the thickness of the 
region through which chemical species must be transported 
in order to be destroyed by nuclear processes. The formation 
of a solar-like convective envelope in a young, differentially 
rotating, 2 Mq star could excite global oscillations, and be 
accompanied by the operation of a solar-like hydromagnetic 
dynamo. Dynamo-generated fields that diffuse into and are 
retained by the radiative interior (see, e.g., Dikpati, Giknan, 
& MacGregor 2006) could enable the star to remain magnetic 
long after spin-down and the elimination of nonuniform rota- 
tion have led to the disappearance of both the surface convec- 



tive layer and the dynamo. Each of these possibilities will be 
addressed in forthcoming papers. 
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